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Recent results from direct detection experiments (Dama, CoGeNT), though subject to debate, 
seem to point toward a low mass (few GeV) dark matter (DM) particle. However, low mass 
DM candidates are not easily achieved in the MSSM nor NMSSM. As shown by some authors, 
singlet extensions of the MSSM can lead to GeV mass neutralinos and satisfy relic abundance 
constraints. We propose here to extract indirect detection constraints on these models in a 
generic way from cosmic-ray anti-proton measurements (PAMELA data) 

Recent results from dark matter direct detection experiments DAMA^ and CoGeNT^l have 
risen a lot of interest for light dark matter (masses between 3 an 20 GeV, roughly). Though 
these results are in tension with those of the XenonlOO^ experiment (in the case of a pure 
spin-independent scattering), the light dark matter hypothesis is an interesting one to explore. 

As it has been shown by Julien Lavalle^, interpreting the results from DAMA^ and CoGeNT- 
in terms of dark matter may be in tension with anti-proton cosmic rays constraints. Indeed a 
dark matter candidate that has such a large coupling to quarks (as suggested by the CoGeNT 
results) may imply a strong annihilation into quark pairs which, may produce more cosmic ray 
anti-protons in the Galactic halo than what has been observed by PAMELA It has been 
shown^that some part of phase of the MSMM may survive the cosmic-ray constraints, however 
only if the astrophysical parameters are very favourable. 

Following the idea of other workslZl§l 9 l 10 l, we have investigated the possibility that rather than 
having a dark matter particle annihilating into lepton pairs, it would annihilate into scalar and 
pseudo-scalar particles which, in turn would decay into Standard Model fermions. There still 
would be a cosmic anti-proton production, however, which a much less sharp spectrum, and, in 
some cases, compatible with the PAMELA measurements. The Next-to-Minimal SuperSymetric 
Model (NMSSM) may have such a particle content (in the case where the dark matter particle 
is mainly singlino it can be made very light and annihilate into the scalar and the pseudo-scalar 
Higgses). However NMSSM is not a requirement and any particle physics model with a similar 
particle content would give comparable results. 

1 Particle content and cross sections 

The phenomenology of singlino-like dark matter scenarios is mostly set by the couplings between 
the neutralino, a Majorana fermion, and the light scalar and pseudo-scalar particles of the Higgs 
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sector. Denoting x the neutralino field and fa the scalar or pseudo-scalar fields, the effective 
Lagrangian that we consider reads: 

£ c ff = - g ^2 x C w x & ~ o ^ ^' ^ k ' 

where we Aj j & is a dimensional coupling taken either fully real or fully imaginary, and where 

^"X* = C XX* + ^XX*T5 
^X* = C XX* ~~ C xx i75 

features both the couplings of the scalar and pseudo-scalar fields fa to the dark matter fermionic 
field x- With this parametrisation, considering the diagrams of Figure [TJ the annihilation cross- 
section can be computed analytically (see details hP}. 

All we have to do is hence to restrict the parameter space of masses and coupling that is 
compatible with the CoGeNT region and collider constraints. 
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Figure 1: Dark matter annihilation channels considered. 



2 Relic Density 

In the standard A-Cold Dark Matter model (ACDM), the dark matter we consider here is ther- 
mally produced during the first times of the universe. Initially in chemical and thermal equi- 
librium with the plasma, its abundance obeys standard quantum statistical equilibrium. When 
the temperature drops below their mass, dark matter particles still annihilate and experience 
Boltzmann suppression until the expansion rate of the universe, controlled by the effective de- 
grees of freedom g+, gets larger than the annihilation rate (see e.g.^). At that time (T ~ ^jy), 
dark matter freezes out, and its relic comoving density is fixed. The present abundance is 
£lh 2 oc -^jy — • It is important to note that when considering low mass dark matter, the de- 
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coupling from the thermal bath occurs at a temperature close to the one of the QCD phase 
transition which strongly affects the value of g±. 

As one can see from Figure [2j for a dark matter particle of 10 GeV, the relic density may 
vary by 60% depending on the QCD phase transition, which is much more than the observational 
uncertainty. This consideration is not new at all, however, as it is more usual to consider higher 
dark matter masses, it is worth recalling. 

Because of the high accuracy of WMAP7 data^, the relic density is extremely constraining 
and as one can see from Figure |4j it reduces considerably the number of possible values of the 
couplings of our model. With this reduced parameter space it is now possible to test the cosmic 
anti-proton flux. 

3 The pbar spectrum 

Before estimating the anti-proton flux at the Earth, one needs to know the anti-proton spectrum 
before propagation, that is the one after the dark matter annihilation into scalar h and pseudo- 
scalar a, their subsequent decay into quark pairs and their hadronization. The first step is easy: 
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Figure 2: Dark matter relic density as a function of dark matter particle mass. The green band corresponds 
to WM AP7 data. The red and blue lines are the results when considering a first order QCD phase transition 

happening at 400 MeV and 150 MeV respectively. 



it is two body annihilation so all the 4-momenta are set by kinematics. In the approximation 
that the dark matter particles annihilate at rest in the halo frame the energy of particle 1 (either 
h or a), is 



4m^. + m\ — m\ 
Am, 



'■x 

and the norm of its momentum is : 



J A (4m| , m\,m^) 



4m x 

which are enough to change from the halo frame to the rest frame of particle 1: (7, f3)={Ei/m\, k/E\) 
In the rest frame of particle 1, the quarks it decays into have energy E* = mi/2 and momentum 

I A;* I = J mf/A — rrig. One finally gets the energy of the quarks in the halo frame: 

E l /m V A (4m|,mj,ml) / A m \ E x 

E q = — - cos(# -1 \ 1 = — - cos(t9)5. 

2 8m x y mi 2 

So the energies of the quarks and anti-quarks coming from the decay of particle 1 are evenly 
distributed between — £ and Finally one gets the probability of having an anti-proton 

of energy Ep from a quark of energy E q f(E q ,E f ) thanks to the P YTHIA0 packaged So finally 
the anti-proton spectrum after the dark matter annihilation is: 

i=l,2 q ^ 2~ £ 

where the first sum is done over the annihilation products (l,2)=(/i, a), and the second sum is 
over all the quark flavours for which 2m q ^ rrtj is satisfied. The branching ratios BRi q depend 
on the particle physics model considered. 



b For this work we made use of version 6.4.24 with CDF tune A 
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Figure 3: Anti-proton flux at the Earth as a function of particle kinetic energy. Data are from PAMELA. Green 
line: background estimation. Brown, yellow and blue correspond to direct annihilation into bb, annihilation into 
a scalar and a pseudo-scalar of masses of 7 and 8 GeV or 3 and 12 GeV, respectively. Dashed are signal only, 

lines are signal plus background. 

Finally, one has to propagate the anti-proton from the annihilation place to the Earth. 
The model used here to describe this propagation has been detailed at length in many other 
papers E 5 I 16 I 17 [ In this model, the charged cosmic rays diffuse off the inhomogeneities of the 
Galactic magnetic field, they interact with the interstellar gas when they cross the disk and 
finally reach the Earth. This model has been shown to be extremely accurate in describing 
many cosmic ray species and, in particular, describes very well the anti-proton astrophysical 
background. This background is due to the spallation of cosmic ray protons and a on the 
interstellar hydrogen and helium (secondary cosmic rays). Moreover, when considering anti- 
protons it is important to also take into account tertiary cosmic rays from inelastic scattering 
of cosmic ray anti-protons. The prediction is in very good agreement with present data (as one 
can see from Figure [3]) and suffers very little from the uncertainties affecting the propagation 
parameters. 

In order to determine whether or not a point in our parameter space is in agreement with 
anti-proton constraints, we simply summed the astrophysical prediction and the dark matter 
component and checked if the total was flux was higher than the PAMELA data. 

4 Results and conclusions 

As one can see from Figure |4j the parameter space which gives a signal in the CoGeNT region, 
gives the correct relic density and in the same time does not give a too high cosmic ray anti- 
proton flux is quite small. This is not the place for a thorough study of the parameter space 
(see more details in^) , however some features can be stressed thanks to Figure |4j 

Large masses of the scalar particle h are forbidden by cosmic ray constraints. Indeed, when 
rrih is too large, the quark spectrum is not very much boosted (the quark pairs are almost 
produced at rest) so the anti-proton spectrum is peaked and easily exceed the observation, as 
if the dark matter particles were annihilating directly into quark pairs. Conversely, very low 
scalar masses cannot be constrained as they cannot decay into bb pairs. From the right panel 
of Figure |4| it appears that, quite naturally, high annihilation cross sections at zero velocity 
(< avo > > 1CP 27 cm 3 .s _1 ) are excluded but not very low ones. The large discrepancy that 
can occur between annihilation cross sections at rest and at time of decoupling shows that the 
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Figure 4: Probing the parameter space. Left: scalar and singlino masses. Right : annihilation cross-section at 
time of decoupling and in the halo. Blue, points in the CoGeNT region, orange: correct relic density but excessive 
anti-proton production, brown, correct relic density and correct anti-proton flux. 

relic density can sometimes been set by the i-channel annihilation rather than by the s-channel, 
alleviating the constraints from anti-proton data. 

Of course, one needs to repeat this work for the different possible dark matter halo profiles 
(here we made use of the profile proposed by^) and for different propagation parameter sets. 
However one can already conclude that interpreting the recent direct detection experiments 
results in term of dark matter is quite challenging. The compatible parameter space, even 
when the dark matter annihilation does not directly goes into quark pairs, is extremely reduced. 
If the CoGeNT, DAMA/Libra result were to be confirmed by CDMS, Xenon and Edelweiss, 
it is interesting to stress that the absence of signal in the cosmic anti-proton channel is very 
enlightening from the point of view of the nature of dark matter as it would put very strict 
constraints on masses and couplings which could be challenging for LHC. 
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